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Very high energy (VHE >100 GeV) gamma rays coming from blazars can produce pairs when 
interacting with the Extragalactic Background Light (EBL) and the Cosmic Microwave Back¬ 
ground, generating an electromagnetic cascade. Depending on the Intergalactic Magnetic Eield 
(IGME) intensity, this cascade may result in an extended isotropic emission of photons around 
the source {halo), or in a broadening of the emission beam. The detection of these effects might 
lead to important constrains both on the IGME intensity and the EBL density, quantities of great 
relevance in cosmological models. Using a Monte Carlo program, we simulate electromag¬ 
netic cascades for different values of the IGME intensities and coming from a source similar 
to 1ES0229H-200, a blazar with hard intrinsic spectrum at redshift z = 0.14, which is an ideal 
distance for potentially observing the effect. We study the possible response of a generic future 
Cherenkov telescope using a simplified model for the sensitivity, effective area and angular res¬ 
olution. Finally, combining these instrument properties, we calculate the angular distribution of 
photons and develop a method to test the statistical feasibility of detecting the effect in the near 
future. 
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1. Introduction 

The Universe is opaque for gamma rays in the VHE (>100 GeV) range. Photon absorption 
in the intergalactic (IG) photon backgrounds is energy dependent and starts to become substantial 
at TeV energies [1]. In particular, VHE gamma rays from jets of active galactic nuclei (AGN) can 
interact with photons in the IR-UV range present in the Extragalactic Background Eight (EBE) 
and photons from the Cosmic Microwave Background (CMB), producing electron-positron pairs. 
These pairs carry most of the energy from the original photons, and can interact as well with IG 
photons via Inverse Compton, promoting them to energies in the HE (>100 MeV)-VHE range, 
making them able to pair produce in the IG backgrounds again. This cascade process converts the 
initial VHE photons into photons of lower energy which can travel further. Moreover, depending 
on the intensity (B) of the Intergalactic Magnetic Eield (IGME), the bending effect on the electron- 
positron pair trajectories can result into different emission scenarios. Eor a strong IGME intensity 
(B > 10^^ G) synchrotron cooling becomes dominant and no secondary gamma rays are produced 
[2]. Eor a moderate IGME (10^^^ G < B < 10^^ G) the electron and positron pair trajectories are 
isotropised around the source eventually giving rise to an extended isotropic emission of photons 
around the source, or halo, which take much longer to reach the observer than the direct photons 
from the source [3]. Eor a weak IGME (B < 10^^^ G) the cascade develops almost exclusively in 
the forward direction, although there is a broadening of the original emission beam, even for very 
small IGME intensities. The extension of this emission depends on the IGME intensity and should 
be clearly distinguished from the halo emission because in this case the broadening takes place 
along the jet direction, not in an isotropic way [4]. Since this effect was proposed several groups 
have tried unsuccessfully to observe it using multiple methods: [5], [6], [7], [8] and [4]. All this 
studies were done using blazars, a subtype of AGN that have their jets pointing towards the Earth 
and are therefore extremely luminous objects, perfect candidates to perform this type of studies. 
The detection of the halo, or the magnetic broadening might lead to important constrains on both 
the IGME intensity and the EBE density, quantities that currently remain uncertain and are of great 
importance in cosmological models. 

In this work we use MC based simulations of intergalactic cascades under different IGME scenar¬ 
ios, and study the effects of the magnetic field on the resulting spectral and angular distributions 
of the aniving photons. Motivated by the existence of a real future Cherenkov telescope system, 
the Cherenkov Telescope Array [9], we assume a simplified model of response for such a type of 
system to develop a method for testing the feasibility of detection of extended emission around a 
given source. 

2. Simulation 

We use the EEMAG program [10] to simulate the intergalactic electromagnetic cascades. This 
simulation takes into account pair production, inverse Compton, synchr'otron losses, and deflections 
of the charged particles in the IGME. We consider a source similar to lES0229-1-200, an HBE type 
blazar at redshift z = 0.14, that seems to have a particularly hard intrinsic spectrum, an important 
feature for the extended emission effect. We inject photons with energy spectrum following a power 
law with index E = 2/3 [11], and adjust the VHE part of the flux to match current observations of 
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this blazar [12]. Here we do not try to explain the SED of this particular source, e.g. we are not 
contemplating Fermi-EAT observations at lower energies, neither considering other values for the 
spectral index. The mention of the blazar 1ES0229-I-200 is just to justify the existence of such a 
spectrum. 

We study the cascade considering different IGMF intensities: B = 10^*^ G, B = 10^^"^ G, 
B = 10^*^ G and B = 10^*^ G, which cover the range of interest for this effect [4]. We also 
consider the emission to be collimated in a cone of 6° [13], which is a more realistic treatment of 
blazar emission than previous approaches. Using this simulations we calculate the SED and the 
angular distribution (0^) of photons arriving at Earth in the energy range 30 GeV - 10 TeV, which 
is the approximate energy window of future Cherenkov telescopes (see Figure 1). 



Figure 1: Spectral energy distribution of photons arriving at Earth as generated by the simulation and 
for different IGMF intensities. The SEDs are scaled to match the VHE part of the measured flux from 
1ES0229H-200 (black points). 


It is interesting to note the differences between the SEDs shown in Figure 1, particularly 
for energies close to 30 GeV (dashed line), where the number of photons increases as the IGMF 
intensity weakens. This is a direct result of the deflection of pairs that end up scattering photons 
in directions off the source-earth line. Differences in the SED can be used to set limits to the 
IGMF ([11], [16]). For an intensity of B = 10^^^ G or weaker, radiation is collimated enough 
to make SEDs indistinguishable. Angular distributions however, present a broadening effect even 
at this weak field regime. Our study focuses on the possible broadening effect that a relatively 
weak IGMF may have in the 6^ distributions, and how this effect can be used to set limits to the 
IGMF when the SEDs show no significant differences. Following this idea, and considering that 
for weaker magnetic fields SEDs will present minimal differences, we continue our analysis for the 
case of B = 10^*^ G. 
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3. Telescope Response 

As mentioned before, there have been several attempts to observe the effect, in all cases only 
upper limits could be obtained. Detection was not possible because the effect is either non existent, 
or too weak to be detected with present instruments. Future generation Cherenkov telescopes will 
enhance their sensitivity and angular resolution improving chances of detection. In this study we 
adopt the generic performance curves proposed by [14]: 

5 = -13.1-0.33X + 0.72x2 (3.1) 

A = 6 + 0.46 A-0.56X2 ^.3 2 ) 

V /68 = 0.038 + exp[-(X + 2.9)/0.61] (3.3) 


where X = log (Photon Energy/TeV), S = Differential Sensitivity/(erg cm^^ s^*), A = Effective 
Area/m^ and t/rgs is the 68 % containment radius of the point spread function (PSE) in degrees. We 
also use the PSE shape obtained by H.E.S.S. [15], 


PSF = C 


exp 



+ Crel exp 



(3.4) 


where Crei = 0.15 and C is a global normalization factor. We combine this PSE shape with t/rgg 
from equation (3.3) to construct the PSE corresponding to a future Cherenkov telescope. The 
resulting PSE is energy dependent and, as expected, narrower than the original in all the energy 
range considered in this study. Eor an event energy of 1 TeV the new function is ~4 times narrower 
than that corresponding to H.E.S.S. (Eigure 2). 



Figure 2: PSF for H.E.S.S. (red) and the estimated PSF for a future array (blue) for Ey=l TeV. 
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A narrower PSF means better angular resolution, which increases the telescope capability of 
distinguishing an extended component. 


4. Method and results 


There are several possible ways for searching extended emission, most of them involve com¬ 
paring the source angular distribution (0^) with the one corresponding to a point source, i.e. the 
telescope PSF. Any significant difference between these two distributions would indicate the exis¬ 
tence of an extended contribution. The ELMAG code generates the 6^ distribution of gamma rays 
when they reach the Earth, however, in order to compare it with the telescope PSE we need the 
6^ distribution as seen by the telescope system. To do this we developed an algorithm that cal¬ 
culates the angular distribution of the averaged number of gamma rays per unit of 6^ detected by 
the telescopes for a given luminosity (L) and observation time {Tobs)- This is achieved by 

sampling the effective area and the PSE for the future array calculated in section 3 and convolving 
it with the 0^ output of the simulation. 

A Cherenkov telescope not only detects events from gamma rays but also from cosmic rays, 
which are the most important source of contamination, or background. In order to discriminate 
gamma-ray from cosmic-ray detections the background must be well understood and properly ex¬ 
tracted from the overall signal. In this case the background rate {b) per unit of 0^ and energy is 
estimated by using the signal {Ny) to noise (Nback) ratio equation, 

Ny ^ (^{E)A{E)AETobs ^4 

^b{E) yI^{E) AE Tobs’ 


where E is the energy, ^(E) is the gamma-ray flux, A(E) is the effective area, and AE is the energy 
bin width. The ability to discriminate background events together with the technical requirements 
to establish a source detection determine the sensitivity of a telescope system. By definition, if 
Ny/s/Nback = 5, then the gamma-ray flux (p (E) corresponds to the sensitivity S{E)/E^. Therefore, 
we have: 


HE) 


'S{E)/E^A{E)H 
5 WesiE) 


AE Tobs- 


(4.2) 


The total background rate {bj) is obtained combining Eqs. (3.1), (3.2), (3.3) and (4.2), and inte¬ 
grating in the considered energy range from 30 GeV to 10 TeV. 

Having the distribution of both the gamma rays from the source and the background rate, the 
next step is to construct the distributions corresponding to the ON and OEE source observations that 
allow to compute the significance of the simulated observation. To construct the ON distribution 
the number of events in the i-th 0^ bin is obtained by taking a random value from a Poisson 
distribution with mean /i (0^)^- = x A0^, where df is the center of the bin and A0^ is 

the bin width. The background in each 0^ bin is calculated also by taking a random value from a 
Poisson distribution with mean jj. (b) =bT x A0^. Then, the ON distribution is obtained by adding 
the background to each 0^ bin. In a real observation, we consider that ON and OFF are taken 
from the same data set, i.e. wobble mode, in which the source is place off the center of the camera 
allowing to have more than one region in the sky where to compute the OFF [17]. We then construct 
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a wobble background considering 3 such regions in the sky (OFF) by taking a random value from a 
Poisson distribution with mean n (b) = 3 x br x Once obtained, this background is extracted 
from the ON distribution to get the excess distribution: ON-OFF/3. 

The comparison between the point source distribution (PSF) with the excess distribution was 
done using the test as follows: 


^2 _ {Nex,i NpsF,i 


i Nex,i + %NoFF,i 


(4.3) 


where Nex.i is the excess and NpsF,i is the weighted sampled PSF in the /-th bin. The denominator 
contains the fluctuations of the quantities involved in the calculation, which follow Poisson statis¬ 
tics. NoFF,i is the fluctuation that comes from adding and subtracting the background (wobble) to 
obtain the excess. Note that the test is applied after both distributions are scaled so that the total 
number of events match to each other (^i^PSFj = 'Li^exj)- 

The test assumes Gaussian statistics for its variables, so a relatively large number of counts 
per bin is required. A common criteria is that up to 20% of bins with < 5 counts is considered 
acceptable for the test to be used [18]. Given the nature of the distributions being compared, where 
most bins contents fluctuate around 0, the validity of the x^ test is not obvious. We therefore cal¬ 
culated the percentage of bins with less than 5 counts which in our case turn out to be significantly 
lower than the requirement (i.e. < 1%). We also calculated the distribution of x'^ and compare it 
with the c/i/-square distribution for a fixed number of d.o.f These distributions show no significant 
differences, suggesting that the c/i/-square distribution is a good approximation of our probability 
function. 

We take the non existence of an IGMF as the null hypothesis. For a fixed source luminosity and 
observation time, our code generates 1000 samples of the excess distribution and uses the test to 
reject the null hypothesis for a fraction of them. This fraction depends on the rejection probability, 
which in this work is taken as 99 %. Figure 3 shows the fraction of samples for which the null 
hypothesis is rejected as a function of the observation time To^,. Although these results appear to 
be optimistic, one has to take into account that the models used to describe the telescope response 
are simple approximations that may be overestimating the capabilities of the telescopes. 


Conclusions 

The IGMF has a direct influence on the development of cascades from distant gamma-ray 
sources, which can be seen in both the spectral energy distribution and the angular (0^) distribution. 
Unlike other previous studies that rely on differences in the SEDs to limit the IGMF intensity, 
our method focuses in the broadening effect that the magnetic field may produce in the angular 
distributions. In this way, the method allows to explore the chance of detecting extended emission 
for weaker IGMF intensities than other methods. Combined with the ELMAG simulation, the 
method can be applied virtually to any VHE source in different redshifts and IGME scenarios, 
offering a way to probe what the response of next generation Cherenkov telescopes would be, and 
what the chances are to detect these effects in the near future. 
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time. 
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